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Abstract 

We construct a deterministic approximation algorithm for computing a permanent of a 0, 1 n by 
n matrix to within a multiplicative factor (1 + e)", for arbitrary e > 0. When the graph underlying 
the matrix is a constant degree expander our algorithm runs in polynomial time (PTAS). In the 
general case the running time of the algorithm is exp(0(r7,3 log^ n)). For the class of graphs which 
are constant degree expanders the first result is an improvement over the best known approximation 
factor e" obtained in |LSW00j . 

Our results use a recently developed deterministic approximation algorithm for counting partial 
matchings of a graph jEGK"*"] and Jerrum-Vazirani expander decomposition method of |JV96j . 

1 Introduction 

A permanent of an n by n matrix A = {uij) is Perm(G) = X^^- ni<i<n '^^^^^ runs over the 

elements of the permutation group on the set 1,2, ... ,n. When A is a zero-one matrix, Perm(G) counts 
the number of perfect matching in the graph corresponding to the adjacency matrix A. Permanent 
is naturally related to the determinant of A (signs in front of products are removed). Yet, while the 
determinant of a matrix can be computed in polynomial time, the problem of computing a permanent 
belongs to the class even when j4 is a zero-one matrix |Val79j . Thus, modulo a basic complexity 
theoretic conjecture, no polynomial time algorithm exists for computing permanents. 

A lot of research was devoted to constructing approximation algorithms for computing a permanent 
of a matrix. The first breakthrough came from Jerrum and Sinclair [JS89j who constructed a randomized 
fully polynomial time approximation scheme (FPTAS) for 0, 1 matrices satisfying certain "polynomial 
growth" condition. This condition relates the ratio of perfect to near-perfect matchings and requires 
that this ratio is at most a polynomial in n. The algorithm is based on a rapidly mixing Markov chain 
which runs on perfect and near perfect matchings. Using this algorithm as a black-box Jerrum and 
Vazirani constructed a randomized approximation algorithm for an arbitrary zero-one matrix, but with 
mildly exponential running time exp(0(n2 log^n)). A recent dramatic improvement was obtained by 
Jerrum, Sinclair and Vigoda |JSV04j . who constructed FPTAS for an arbitrary matrix with non-negative 
entries. 

Unfortunately, the randomization aspect of the algorithm of [JSV04j is quite crucial. It is not known 
how to derandomize their algorithm and the best known deterministic approximation algorithm was 
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given by Linial, Samoridnitsky and Wigderson [LSWOOj . who provide only multiplicative approxi- 
mation factor. Their algorithm is based on a FPTAS for a related matrix scaling problem and uses 
van der Waerden's conjecture, which states that a permanent of every doubly stochastic matrix is at 
least n!/n". The approximation factor can be improved to {k/{k — 1))^^ for the case of of matrices 
with row and column sums bounded by k, using the Schrijver's bound |Sch98] . For generalizations and 
simplified proof of Schrijver's bound see Gurvits |Gur06j . Thus approximation of a permanent is one 
of the famous algorithmic problems where the existing gap between the randomized and deterministic 
algorithms is so profound. 

In this paper we establish the following two results. First we construct a polynomial time algorithm 
which for every e provides (1 + e)" multiplicative approximation factor for a permanent of a 0, 1 matrix, 
when the underlying graph is a constant degree expander. The definition of the latter is given in 
the subsequent section. While our algorithm is polynomial, it is not fully polynomial, as the term 1/e 
appears in the exponent of the running time. Thus we significantly improve the e"' factor of [LSWOOj for 
this class of graphs. Next we construct an algorithm providing the same approximation factor (1 + e)"" 
for a permanent of an arbitrary 0, 1 matrix. The running time of the algorithm exp(0(n3 log'^ n)). The 
main technical ingredient of our results is the reduction of the problem to computing a partition function 
(see the next section for the definition) corresponding to the collection of all partial and full matchings 
of the graph. Recently there was a progress in constructing deterministic approximation algorithms 
for various counting problems without the usage of Markov chain based methods. The approach was 
introduced by Bandyopadhyay and Gamarnik [BG06j . and Weitz [WeiOGj and is based on correlation 
decay property originating in statistical physics literature. Subsequent works in this direction include 
Gamarnik and Katz |GK07j and Bayati et al. [BGK"*"] . It is the second of these two works which 
provides the basis of our work. Specifically, a polynomial time approximation algorithm is constructed 
in BGK"*" for computing the partition function corresponding to partial matchings of a an arbitrary 



constant degree graph. The same algorithm solves the underlying problem for an arbitrary graph but 
with the running time exp(0(n2 log^n)). In this paper obtain appropriate bounds on the permanent 
of a matrix in terms of the partition function of partial matchings of a graph. Additionally, our 
second result corresponding to arbitrary graphs relies heavily on the expander decomposition technique 
developed by Jerrum and Vazirani |JV96j for constructing a mildly exponential approximation algorithm 
for computing a permanent of a 0, 1 matrix. 

The remainder of the paper is organized as follows. Definitions, preliminaries and the statements of 
the main results are subject of the following section. Section [3] is devoted to constructing approximation 
algorithm for the case when the underlying graph is a constant degree expander. The general case is 
considered in Section SI We conclude with some possible further directions in Section [5l 



2 Preliminaries and the main result 

Consider a simple undirected n by n bi-partite graph G with the node set V = ViL)V2, \Vi\ = | V2I = n. 
Let E be the set of edges of the graph. N{v, G) C ^ denotes the set of neighbors oi v £ V and 
A(f) = |A^(t',G)| denotes the degree of the node v. The degree of the graph is A = max^, A(t>). Given 
a set of nodes A, we denote by N{A) or specifically by A'^(^, G) the set of neighbors of A (in G). Given 
a > 0, a graph G is defined to be an expander or specifically a- expander if for every set of nodes 
A cVi,i = 1,2 such that |^| < n/2, the inequality \N{A)\ > (1 + a)\A\ holds. We also define 

a{G) = max — — 1 

A A 
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to be the expansion of the graph G, where the maximum is over all subsets A C Vi,i = 1,2 with 
1^1 < n/2. Clearly, G is a-expander if its expansion is at least a. 

A matching is a subset M <Z E such that no two edges in M share a node. For every k < n let 
M{k) be the number of size k matchings in G. Specifically, M{n) is the number of full matchings. 

Given a graph G\et A = (uij) be the corresponding adjacency matrix. The rows and columns of A 
are indexed by nodes of Vi and V2 respectively, and aij = 1 if {vi,Vj) G E and aij = otherwise. It is 
immediate that M{n) = Perm(A). 

A parameter A > 1 is fixed called the activity. The partition function corresponding to A is defined 

as 

Z(A,G) = J^AI^^' = Yl ^''M{k). 

M l<k<n 

A partition function is an important object in statistical physics. The case of matching is usually called 
monomer-dimer model in the statistical physics literature. 

Definition 1. An approximation algorithm A is defined to be a Fully Polynomial Time Approximation 
Scheme (FPTAS) for a computing Z(A,G) if given arbitrary 5 > it produces a value Z satisfying 

1 Z 
TTs - Z{X,G) - ^ + 

in time which is polynomial in n and ^ . 

The following result was estabhshed in [BGK+j . 

Theorem 1. There exist a deterministic algorithm which provides a FPTAS for computing Z{X,G) for 
an arbitrary graph/activity pair G, A when A and A are constants. The complexity of the same algorithm 
is exp(0(\/Anlog^ n)) for general A and A. 

The case A = 1 corresponds to counting the number of partial matchings of a graph. In this paper 
we use the algorithm underlying Theorem [T] as a subroutine to devise an approximation algorithm for 
computing a permanent. For this purpose we will be making A to be appropriately large. Throughout 
the paper we assume A > 10. 

We now state the main two results of this paper. In the next and the following results the notion of 
(1 + €)"■ multiplicative approximation factor of Perm(G) corresponds to obtaining a value Z satisfying 

Theorem 2. Let G be an n by n bi-partite a-expander and let e > 0. There exist a deterministic 
(1 + e)" approximation algorithm for computing Perm(G) with running time exp{0{V e^^a^^n log^ n)). 
Moreover, the running time is polynomial in n whenever A, a are constants. 

Our algorithm corresponding to the second part of the theorem, while polynomial, is not strongly 
polynomial. As we shall see, the dependence of the running time on the approximation parameter e 
is of the form n^^~\ While our approximation factor (1 + e)" is a far cry from PTAS (namely 1 + e 
approximation factor), it is still a significant improvement over e" factor constructed in [LSWOOj . 

Our second result does not require any restrictions on the underlying graph. 

Theorem 3. There exist a deterministic (1 + e)*^ approximation algorithm for computing Perm(G) of 

12 o 

an arbitrary n by n bi-partite graph G which runs in time exp(0(e~2n3 log n)). 

Thus, similarly to | JV96j . our algorithm provides a mildly exponential algorithm for approximating 
a permanent (with a weaker approximation factor (1 + e)"). 
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3 Constant degree expanders 



Proof of Theorems [2] is given in this section. We begin by estabhshing some prehminary results. We 
assume without the loss of generality that M{n) > 1. Consider an arbitrary k matching M between 
sets Ai cVi,A2 C V2, l^il = 1^2! = k. A path vi,V2, ■ ■ ■ , V2r is defined to be an alternating path wrt M 
if vieVi\ Ai,V2r G ^2 \ A2, if iv2i,V2i+i) G M, 1 < Z < r - 1 and {v2i-i,V2i) ^ M, 1 < / < r. Observe 
that given M and an alternating path P one can construct a k + 1 matching, by subtracting from M 
all the edges {v2i,V2i-{-i), 1 < I < r and adding all the edges (f2Z-i, '^2;)) 1 < ^ < r — 1. The length of this 
alternating path is defined to be 2r. 

Lemma 1. Let AI be a k matching between Ai and A2 for k < n — I. For every set L C Vi\Ai with 
\L\ > {n — k)/2 there exists an alternating path P with end points in L and V2 \ A2 and with length at 
most 

o{^og{-^)\og-\l + a)). (1) 

Proof. The proof is similar to the argument of Lemma 2 |JV96j . Let R = V2 \ ^2- Let Lr be the 
set of nodes in V2 reachable from L en route of alternating paths with length at most 2r. Let tq 
be a minimum r satisfying (1 + a)^(n — k)/2 > n. Note = log(^^^) log^^(l + a)j . Since 

\L\ > (n — k)/2, then /(rg) — min((l + a)''°|L|, ^ + 1) = ^ + 1. For every r < tq, by the expansion 
property either \Lr\ > min((l + a)^|L|, ^ + 1), or L^i n i? 7^ for some r' < r. In the second case we 
found an alternating path with length < ^(ro). In the first case we have |-Lrol > n/2. We now claim 
that for every r > ro either n 1^2 \ ^2 / or |y2 \ -^^r| < -jj^^^yF^in/^). Indeed, let L'^. C Vi be the 
set of nodes matched to Lr- In particular \Lr'\ > n/2. Since |-^^r+i| = \^{Lr')\ > n/2, then 

\V2\Lr\ = \Vi\Lr'\ > \N{V2\Lr+l)\ > {l+a)\V2\Lr+l\, 

and the assertion is established by induction. For r > 2ro we obtain ^-^_^^Y-rQ {n/2) < n — k and thus 
we must have H V2 \ A2 7^ 0. We established that for some r < 2ro there exists an alternating path 
with an end points in L and V2 \ A2. The length of this path is 0(ro) = 0(^^log~^(l + a)). This 
completes the proof. □ 



Lemma 2. For every k < n — 1 

M{k + 1) ~ ~Vn 



M{k) / n No(^iog-i(i+«)iogA) 

Uk + l) - \n-k) 



(2) 



As a result 



M{k) ^ / 2en \ o(^(n-fc)log-i(l+a)logA 
M{n) ~ \n - k) 



(3) 



Proof. Fix an arbitrary k matching M between Ai C Vi and A2 <Z V2. We claim that there exist at 
least {n — k)/2 k + 1-matchings obtained from M via an alternating path with length at most the value 
given by ([T]). Indeed, consider the set of all nodes v \nVi \ Ai such that the shortest alternating path 
starting from v is larger than the required bound. By Lemma [T] this set contains less than (n — k)/2 
nodes and the assertion follows. 

Now consider the following bi-partite graph. The nodes on the left (right) are all k (A;+l)-matchings. 
We put an edge between two matchings M, M' if M' is obtained from M via an alternating path with 
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length at most ([T]). The total number of edges in this graphs is at least M{k){n — k)/2 by our observation 
above. For every matching M' on the right side of the graph the number of edges pointing to it is at 
most the number of alternating paths with length at most ([1]) which result in M'. For every possible 
starting node of an alternating path, the number of such alternating paths is at most 



^0(^log(^)log-i(l+a) 

The number of starting nodes is bounded by n. Then the total number of edges in this graph is at most 

o( log(^)log-i(l+a) ) 

M(fc + l)nA V /. We conclude that 

of log(^)log-ni+a) I 

M{k + l)nA V ^^"-^^ ^ ^ 'J >M{k){n-k)/2. 
The bounds ([2]) then follows. From this bound we also obtain 

M{k) <2"-fcf """'^ _/ (2en)"-'= x o(^iog-i(i+Q)iogA) 

M(n) - \{n-ky.) ~ \{n-k)^-'') 

where the Stirling's approximation was used in the equality. This is the required bound ([3]). □ 
Corollary 1. The following holds 

1 < ;,n^ptm(G) ^ exp(0(nA-i log-i(l + a) log A)). (4) 

Proof. The inequality Z{X,G) > A"Perm(G) is immediate. We focus on the second inequality in (jj]). 
We need to analyze the ratio 



X'Mjk) _^_(^_,)M{k) 
\^M{n) M{n)' 



(5) 



We set Cn = log (1 + a) log A for simplicity. Applying the second part of Lemma [2] 

A-("-')^ - ^""P {0{cn{n -k){\ogn- log(n - k) - log ^ 

Consider the problem of maximizing 

g{x) = X log n — X log x — x log — 

2e 

in the range x G [l,n]. The boundary cases x = l,n give respectively values log(2en/A), — nlog(A/2e). 
The second quantity is negative when A > 2e (recall our assumption A > 10). To find another candidate 
for the largest value, we take the derivative with respect to x and equating it to zero we obtain 

log n — log X — log = 

giving X = 2n/\. Evaluating g{x) at this value simplifies to 2n/A and this gives the largest value of g 
when n is larger than some A dependent constant. We conclude that the left-hand side of ([5]) is at most 
exp(0(^)), implying 

Z{\,G) f/^r^niT-s. //o/''^""' I 1 ^^ 

< nexp(C'(— — )) = exp(C'(— — + lognj) = exp^O^— — )). 



A"Perm(G) ~ ' A ' X ° " ' X 

□ 
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Proof of Theorem\^ Fix an arbitrary constant e > 0. We select the smallest A so that exp(0(A^^ log^^(l+ 
a) log A)) < 1 + e. It is clear that A = 0{e~^a^^ log A). We compute an e- approximation Z of Z{\, G) 
using an algorithm from Theorem [T] for computing Z'(A, G). By Corollary [U it satisfies 

Then Z/A" provides the required approximation. The complexity of this algorithm is 

exp(0(\/An log^ n)) = exp(0(\/ e^^a^^ (log A)n log^ n)) = exp(0(\/ e^-'^a'-'^nlog^ n)), 

and the first part of the theorem is established. 

For the second part we observe that A = 0(e^^a^^ log A) is a constant whenever a and A are 
constants. We recall from Theorem [T] that the algorithm for computing Z(A,G) is polynomial time, 
under these assumptions. □ 



4 General graphs 

Proof of Theorem O Our approach borrows heavily from the Jerrum-Vazirani expander decomposition 
approach |JV96j . The idea is to decompose the underlying graph into a collection of subgraphs with a 
suitable expansion properties and apply an algorithm for computing a permanent recursively. In [JV96] 
the subroutine used is based on the algorithm relying on rapidly mixing Markov chain. Here we use the 
deterministic algorithm constructed in the proof of Theorem [2j 

The following result is established in [JV96] (Lemma 4). There exists and algorithm, called TestExpansion 
which on input G, a either identifies that G is an a-expander, or identifies a set A C Vi, |^| < n/2 such 
that N{A) < (1 + 2a)|^|. The running time of the algorithm is exp(0(anlogn)). 

Given an arbitrary set A dVi'it is straightforward to observe that 

Perm(G) = ^ Perm(^, B) Perm(A^ B") (6) 

BCN{A),\B\ = \A\ 

where A^ = Vi \ A, B'^ = V2 \ -B and Perm(A, B) is the permanent of the subgraph induced by 
A and B. Based on this observation we propose the following recursive algorithm for computing 
Perm(G). We suppose that we have an algorithm Ar which computes (1 + e)*" factor approxima- 
tion of a permanent of any r by r bipartite graph in time g{r), for every r < n — 1. We use 

it to construct An as follows. Set a = n~^^^. Run algorithm TestExpansion on G. The run- 

2 

ning time is exp(0(n3 logn)). If the algorithm returns no set A, then the underlying graph is an 
a-expander and we use algorithm of Theorem [2] to obtain an (1 + e)" approximation of Perm(G). 

The running time is exp(0( Ve^^nsnlog'^ n)) = exp(0(e~2n3 log'^ n)), and the overall running time 
is exp(0(n3 logn)) + exp(0(e~2n3 log^n)) = exp(0(e~2n3 log'^n)). Let cq be the constant hidden in 
O(-). From now on, treating e as constant, we hide e~2 in the O(-) term. In the end we observe that 
the dependence of the running time on e is of the form exp(0(e~2)). 

Otherwise the algorithm TestExpansion identifies a set A with |A^(yl)| < (1 + 2a)|^|. We estimate 
Vexm^A^B) and Perm(^'^, S'^) using algorithm Ar with r = |j4|,|A'^| respectively. Then we estimate 
Perm(G) using the expression ([6]). For every product Perm(A, i?) Perm(y4'^, S"^) our approximation 
factor is (1 + e)l"^l(l + e)'"^"' = (1 + e)" by the recursive assumption. 
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Now we obtain an upper bound on g{n) and specifically show that it is exp(0(n3 log n)). Let c > cq 
be a very large constant. By the recursive assumption g{r) < exp(cr 3 log^ r)), r < n — 1. The function 
g satisfies the following bound 

^(n) < max f expfcona log^ n), max ( ^ ~'~ ] (g(r) + g(n — r))) . (7) 

V l<r<n/2 V / ^ 

Here the term (J'^^^'^"'^) comes from performing the computation over all subsets B of N{A) with size 
\B\ = \A\. This term is trivially upper bounded by 

^2ar _ gxp(2n^3r logn) (8) 

2 

< exp(n3 logn). (9) 

Now, by the recursive assumption we have for r < n/2 that g{r) < exp(cr » log^ r)) < exp((3/4)cn3 log^ n)), 
which gives 

rfl + 2a)\ 2 

]g(.r) < exp((7/8)cn3 log^n)), 
r J 

provided that 1 + (3/4)c < (7/8)c. As for g{n — r), we have for sufficiently large n 

r(l + 2a)\ / _i , o, 

(7(n — r) < exp(2n 3r log n) exp(c(n — r) 3 log in — rj) 

r J 

< exp(2n~3r logn) exp(c(n — r)3 log n) 

(a) ^ 2 . 1 . 

< exp(2n~ 3 r logn) exp (cn 3 log'^ n — (cr/2)n~ 3 log'^ n) 

. 2, ^ -In 

< exp(cn3 log n — 2n 3 j 

< (1 — n~3 ) exp(cn3 log'^ n), 

2 2 1 2 

where in (a) we use (n — r)3 < ns — |n~3 (obtained for example using Taylor's expansion around ns); 
and in (b) we use r > 1 and c > 8. Using exp((7/8)cn3 log^n)) < n^3 exp(cn3 log^ n) for large n, we 
obtain that for every 1 < r < n — 1 

r(l + 2a)\ , / s , / ^ .23 

' (Sl^j + 9{iT' ~ ^)) ^ exp(cn3 log nj. 

Combining with d?]) we obtain the required bound g{n) < exp(cn3 log^n). □ 



5 Conclusions 

We proposed a new deterministic approximation algorithm for computing a permanent of an n by 
n 0, 1 matrix. Our algorithm provides a multiplicative approximation factor (1 + e)" and runs in 
polynomial time for a matrix corresponding to a constant degree expander, and in time exp(0(n3 log'^)) 
for an arbitrary matrix. Our algorithm is based on a recent deterministic approximation algorithm for 
counting the number of partial matchings of a graph [BGK^j . It is natural to try to extend our results 
in several directions. One possibility is lifting 0, 1 requirement. This entails extending the result of 
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BGK+] to the case of weighted graphs. It is reasonable to expect that for some special cases such a 
program will work. In this case one can try to extend the results of the present paper to the case of 
general matrix entries. 

A second interesting direction is utilizing the ideas in |JSV04j , developed in the Markov chain setting, 
as a possibility of getting stronger deterministic approximation algorithms for computing a permanent. 
The technique developed in [JSV04| allow one to deal with the case when the ratio M{n — 1)/M{n) is 
exponentially large. Since obtaining amenable bounds on the ratios M{k)/M{n) was required for our 
algorithms to work, this might be indeed a fruitful direction for further research. 
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